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I  INTRODUCTION 


There  has  been  substantial  effort  to  model,  and  subsequently  pre¬ 
dict,  the  ionospheric  effects  of  nuclear  explosions  at  high  altitudes. 

The  effects  which  will  have  serious  impact  on  both  satellite  and  ground- 
based  communications  include  sudden  enhancement  of  electron  density,  and 
sudden  heating  of  both  neutral  and  charged  species  resulting  in  atmo¬ 
spheric  heave. 

Large  electric  fields  of  magnetospheric  origin  have  substantial 
influence  on  the  state  of  the  ionosphere.  These  electric  fields,  through 
the  process  of  joule  heating,  produce  dramatic  increases  in  the  ion  and 
neutral  temperatures.  In  some  instances  the  heating  is  sudden,  and 
results  in  an  upwelling  or  heave  of  the  neutral  atmosphere.  This  results 
in  a  change  in  ion  chemistry,  so  that  a  depletion  of  electron  density 
occurs . 

Impulsive  heating  events  have  been  observed  by  the  Chatanika  inco¬ 
herent  scatter  radar  (Chatanika,  Alaska).  Increases  in  ion  temperatures 
of  a  few  thousand  degrees  kelvin  occur.  Whether  the  heat  source  is  a 
distant  nuclear  explosion  or  a  large  joule  heating  event,  the  resulting 
ionospheric  modification  should  be  similar.  However,  we  find  that  the 
nuclear  effects  code,  WEPH  VI  for  example,  docs  not  predict  similar 
effects.  In  fact,  given  sufficient  time  following  a  nuclear  "event"  the 
code  never  returns  all  of  the  parameters  to  ambient  values. 

The  intent  of  the  code  is  to  predict  only  large  scale  perturbations. 
As  a  result,  the  less  dramatic,  but  still  significant,  effects  of  heating 
the  ionosphere  outside  of  the  fireball  arc  generally  ignored.  These 
effects,  though,  affect  a  large  volume  of  ionosphere  changing  the  ion 
chemistry,  and  thus  result  in  sufficient  ionosphere  depletion  to  limit 
HF  propagation. 
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I I  BACKGROUND 


The  high-latitude  thermosphere,  in  contrast  with  the  midlatitude 
thermosphere,  is  strongly  coupled  to  the  magnetosphere.  As  a  result, 
large  amounts  of  energy  and  momentum  are  transferred  from  the  magneto¬ 
sphere  to  the  ionosphere.  These  energy  inputs  have  a  significant  effect 
on  the  morphology,  thermodynamics,  and  chemistry  of  the  high-latitude 
ionosphere.  In  some  cases,  the  magnitude  of  the  energy  deposition  is 
similar  to  that  of  a  distant  nuclear  explosion,  i.e.,  the  temperature 
enhancements  are  similar  to  those  expected  from  nuclear  sources.  Ener¬ 
getic  auroral  particles  are  responsible  for  both  heating  effects  and  the 
production  of  ionization.  Electric  fields  arc  responsible  for  joule 
heating.  The  relationship  between  the  auroral  processes  and  their  mani¬ 
festation  in  the  high-latitude  ionosphere  is  outlined  in  Figure  1. 


FIGURE  1  AURORAL  BLOCK  DIAGRAM  (adapted  trom  Rees,  1975) 
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A .  Particle  Precipitation  and  Electron  Temperature 

On  the  left  side  of  the  auroral  diagram  we  indicate  the  heating 
process  due  to  auroral  electrons.  A  flux  of  energetic  electrons  bombards 
the  upper  atmosphere,  the  depth  of  penetration  depending  on  the  electron 
energy.  Most  of  the  energy  is  deposited  between  100-  and  120-km  altitude. 
The  energetic  primary  electrons  transfer  energy  through  the  ionization  of 
the  ambient  neutral  constituents  (N_ ,  0  ,  and  0).  The  collision  reactions 
arc : 


n2  +  O  - ►  N2  +  2c  (1) 

N2  +  e  - ►  N+  +  N  +  2e  (2) 

02  +  e  - ►  0;  +  2e  (3) 

0  +  c  ►  0+  +  O  +  2e  (4) 

4 

0  +  e  ■■  ■»  0f  +  2c  (5) 


In  the  above  ionizing  reactions  some  of  the  energy  from  primary 
electrons  goes  into  the  production  of  excited  states  of  N+ ,  0+ ,  N,  and  0. 
Enhanced  vibrational  and  rotational  levels  of  N'2  ,  0(?  ,  and  0*  arc  also 
produced.  The  secondary  electrons  have  energies  of  tens  of  electron 
volts.  Much  of  this  energy  is  lost  in  inelastic  collisions  with  neutrals 
and  finally  the  remaining  energy  is  shared  with  the  ambient  electrons 
producing  a  temperature  enhancement. 

15 .  Joule  Heating  and  Ion  Temperature 

Referring  again  to  Figure  1,  the  right  side  indicates  the  process 
of  heating  by  electric  fields  perpendicular  to  the  magnetic  field. 

Magnetospheric-K  fields  are  mapped  down  along  magnetic  field  lines  and 

-¥ 

are  consequently  applied  to  the  ionosphere.  The  field  produces  E  X  B 
drift  of  both  ions  and  electrons  As  the  ions  drift  through  the  neutral 
atmosphere  they  collide  with  neutrals  and  frictional  healing  (or  joule 
heating)  results. 
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The  ion  temperature  can  exceed  both  the  neutral  and  electron  tem¬ 
peratures  at  ionospheric  heights.  Temperature  enhancements  ol  2000  K 
(3  times  ambient)  have  been  observed  as  the  result  of  heating  by  a  60-mV 
orthogonal-E  field. 

C .  Ionospheric  Heating  From  Nuclear  Effects 

The  radiation  from  high  altitude  nuclear  explosions  will  have 
similar  effects  on  the  upper  atmosphere  as  auroral  particle  precipitation. 
The  UV  radiation  from  the  debris  kinetic  energy  and  X-ray  from  prompt  and 
delayed  radiation  arc  the  primary  sources.  The  prompt  and  delayed  radia¬ 
tion  is  less  important  at  F-region  altitudes.  Figure  2  indicates  the 
ionization  process  and  the  resultant  heating. 


▼ 


D .  Auroral  Zone  Ion  Composition 

In  addition  to  the  ionization  reactions,  the  principal  auroral 
ionosphere  reactions  are  as  follows: 


N2  +  ° 


NO  +  N 


N2  +  °2  — *  0+2  +  N2 


N+  +  0„ 


N+  +  0„ 


°2  +  N2 


°2  +  N 


0+  +  N, 


N2  +  O 


°2  +  ° 


NO  *  e 


NO  +  0 


°2  +  N 


NO  +  NO 


N0+  +  0 


0^  +  NO  - ►  N0+  +  02 


N0+  +  N 


0+  +  0„  - ►  0^  +  0 

2  2 


N  +  N 


-*•  0  t-  0 


N  +  0 


(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

(15) 

(16) 
(17) 


■1-2  2 

There  are  also  reactions  involving  excited  states —  0  (  D,  P), 
O^la*  )  which  produce  N*  by 


0+  (2D)  f  N 


0+  (2P)  +  N 


°2  (al  >  f  N2 


N+2  ^  0 


N2  +  ° 


N2  +  °2 ' 


(18) 

(19) 

(20) 
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There  are,  of  course,  other  reactions,  but  those  listed  above  are  the 
principal  ones. 


We  have  investigated  the  effects  of  large  E  fields  on  the  concen¬ 
trations  of  N0+  and  0+  at  high  latitudes.  The  measurements  indicate  that 
during  periods  of  large  E  fields  the  dominant  F-region  ion  becomes  NO+ 
instead  of  0+ .  As  a  result  of  this  change  in  ion  composition,  the  F- 
region  ionization  rapidly  decays  because  of  recombination  of  the  molecu¬ 
lar  ions. 

The  following  processes  contribute  to  increasing  the  N0+  abundance; 
these  are: 

•  enhanced  LN^] 

•  enhanced  ion  temperatures 

•  enhanced  NQ  vibrational  temperatures. 

Enhanced  [N^J  arises  because,  during  periods  of  large  joule  heating, 
it  is  expected  that  the  neutral  atmosphere  should  be  heated  by  ion-neutral 
collisions.  Hays,  ct  al.,  [1973],  have  shown  that  a  heated  volume  element 
of  neutral  particles  will  rise  with  a  vertical  velocity  which  increases 
with  altitude.  This  expansion  will  cause  an  upwclling  of  N2  molecules, 
thereby  increasing  the  number  density  of  at  higher  altitudes.  In  the 
nuclear  case,  this  phenomena  is  called  "heave.”  The  number  density  of 
N0+  will  be  enhanced  because  more  1S'9  is  available  to  cither  react  with  0+ 
to  form  NO^  (Ec; .  13),  or  be  ionized  to  form  (Eq .  1),  which  in  turn 
reacts  with  0  to  form  N0+  (Eq .  6). 

The  second  contributing  factor  to  the  enhancement  of  N0+  ions  is  the 
dependence  of  certain  reaction  rate  coefficients  on  ion  temperature.  The 
rate  coefficient  (K^)  for  the  reaction  0+  +  N2  -»  N0+  +  N  from  Banks  et 
at  .  ,  1 1974 J  ,  is 


K1  =  < 


,  ,  .  n-  12  /  300  K\  3 

1 . 2  X  10  I  — -  J  cm  /  s 

V  eff  / 


8.0  x  10 


-I'M  Teff 


300  K ) 


3, 

cm  /  s 


1  1 


eff 


750 


eff 


..  >  750 


■+/  2  "V  •+ 

where  T  „  =  T  +  0.329  E  ,  and  E  is  the  orthogonal  E-field  measured 
eff  n  x  ± 

in  the  neutral-wind  frame.  This  is  a  simplified  form  assuming  v  A1  « 1 

and  B  =  0.5  gauss  (Polar  region).  As  an  example,  for  an  E  field  of 

50  mV/m,  and  T  =  900  K,  T  will  be  increased  by  a  factor  of  1.9,  as 
n  eff  J  ’ 

compared  to  T  __  for  no  E  field.  This  in  turn  increases  K,  by  a  factor 
eil  1 

of  3.7  .  The  result  should  be  an  enhancement  in  the  number  of  NO  4  ions. 
The  ions  are  lost  through  dissociative  recombination.  The  recombination- 
rate  coefficient  decreases  with  increasing  electron  temperature  [Walls 
and  Dunn,  1974J ;  consequently,  the  net  result  is  an  increase  in  the 
number  density  of  molecular  ions. 

The  third  point  relative  to  the  production  of  NO1  concerns  the 
+  *  +  * 

reaction  0  r  N  *♦  NO  +  N,  where  N  is  vibrational ly  excited  N  .  The 

£  £ 

reaction  rates  for  this  reaction  arc  enhanced  as  compared  to  0  +  N 

[Banks  ct  al . ,  1974] .  In  the  auroral  oval,  N2  becomes  vibrationally 
excited  through  the  process 

* 

c  (energetic  electrons)  f  N  -♦  N  +  c  (21) 

N  +  NO  -»  N*  +  O  (22) 

o  (LD)  +  N2  -»  N*  t-  0  C5P).  (23) 

Alt  of  these  processes — i.e.,  (1)  enhanced  number  density  of  N2 , 

(2)  enhanced  reaction  rate  coefficients  due  to  elevated  T  ,  and  (3) 
enhanced  reaction  coefficients  due  to  vibrational  excitation  of  h'2,  act 
jointly  to  increase  the  N0f  abundance  in  the  F  region.  The  first  and 
second  processes  are  mainly  a  result  of  joule  heating,  and  the  third 
process  results  from  electron  precipitation. 

The  nominal  lifetime  of  molecular  ions  in  the  presence  of  moderate 
electron  densities  (10°  el/cm'*)  is  only  a  few  seconds,  whereas  the  04 
lifetime  is  minutes.  Consequently,  changing  the  dominant  E-region  ion 
to  NO  *  results  in  a  depletion  of  electron  density.  During  the  event 
examined  in  the  study,  the  critical  frequency  decreased  to  about  2.7  MHz 
from  4.7  MHz. 
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Ill  MEASURED  AURORAL  ZONE  IONOSPHERIC  RESPONSE 


The  geomagnetic  coordinates  of  the  Chatanika  incoherent-scatter 
radar  are  given  in  Table  1,  and  its  parameters  are  given  in  Table  2. 


Table  1 

GEOMAGNETIC  COORDINATES  OF  THE  CHATANIKA  INCOHERENT -SCATTER  RADAR 


Geographic 

Coordinates 

Dipole  Geomagnetic 

Coordinates 

Magnetic  Field 

Latitude 

Longitude 

Latitude 

Longitude 

Dip  Angle 

Declination 

65.1 03 °N 

147 .451°\V 

64 . 75°N 

1 05 . 0°\V 

76.5° 

29.0° 

Table  2 

PARAMETERS  OF  THE  CHATANIKA  INCOHERENT -SCATTER  RADAR 


Operating  frequency 

1290  MHz 

Transmitted  peak  power 

3  to  4  MW 

Pu 1 scwidths 

60  t :  s ,  160  us  ,  320  ns 

Effective  antenna  aperture 

180  m2 

Antenna  on-axis  gain 

47.1  dB 

Antenna  1/2 -power- fit  11 -width 
beamwidth 

0.6° 

Transmit  polarization 

Right  circular 

Receive  polarization 

Left  circular 

System  noise  temperature 

110  K 

A/D  converter  sample  spacing 

10  ,;S 

On-line  computer  system 

XDS  930 

•j 
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The  geometrical  relationship  between  Chatanika  and  the  auroral  oval 
for  moderate  geomagnetic  activity  is  shown  in  Figure  3.  The  oval  is 
usually  far  north  of  Chatanika  at  local  noon  and  is  often  overhead  near 
midnight.  The  circle  (arbitrarily  located  at  1500  UT)  indicates  the 
E-region  coverage  typical  for  these  experiments. 

The  Chatanika  radar  measures  the  following  ionospheric  parameters: 

(1)  Electron  density  (N^)  (backscatter  power) 

(2)  Ion  velocity  (VO  (doppler  shift) 

(3)  Temperature  ratio  (T  /T.) 

...  .  (spectral  shape) 

(4)  Ion  temperature  1 

From  these  measurements  other  parameters,  such  as  the  electric  field, 
currents,  and  conductivities,  are  derived. 


oo 

FIGURE  3  CHATANIKA  LOCATION  (65°N)  AND  FELDSTEIN’S 


1967  AURORAL  OVAL.  (Source:  Feldstein  and 
Starkov,  1967.) 
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A. 


Radar  Measurements 


A  representative  joule  heating  event  of  interest  occurred  on  13 
August  1976.  The  radar  was  operated  continuously  for  the  entire  day. 

The  operating  mode  was  an  azimuth  scan  at  a  fixed  elevation  angle  of 
76.5°  (the  magnetic  dip  angle).  The  time  for  one  antenna  scan  was  about 
6  minutes,  and  the  data  were  integrated  over  two  antenna  scans  giving 
12-minute  averages  of  all  the  parameters. 

The  ambient  electron  density  profile  is  shown  in  Figure  4.  This 

profile  was  measured  about  45  minutes  prior  to  a  large  joule  heating 

5  3 

event.  It  can  be  seen  that  an  auroral  E-region  of  about  3  x  10  el/cm 

5  3V 

and  a  moderate  F  region  (2  X  10  el/cm  )  existed. 

The  electric  field  was  computed  for  the  24-hour  period  and  it  shows 

the  characteristic  pattern  of  auroral  zone  electric  fields.  In  the 

evening  prior  to  local  midnight,  the  field  is  directed  to  the  north. 

Near  midnight  the  field  reverses  and  is  directed  south  during  the  early 

morning.  The  north-south  field  is  plotted  in  Figure  5.  The  field 

becomes  quite  largo  at  1400  UT — nearly  60  mV/m.  The  resulting  joule 

heating  is  plotted  in  Figure  6.  At  1400  UT  the  energy  input  is  about 
2 

30  ergs/cm  -s ,  which  is  substantial.  This  is  one  of  the  largest  joule 
heating  events  analyzed  in  detail  and,  as  will  be  seen,  its  effects  are 
considerable . 

The  electron  and  ion  temperatures  at  277  km  altitude  are  shown  in 
Figure  7.  At  the  start  of  the  run,  0000  UT ,  the  ion  temperature  is 
fairly  constant  with  a  value  of  900  K.  There  is  little  change  in  the 
ion  temperature  until  about  1000  UT ,  when  the  joule  input  begins.  There 
is  an  obvious  enhancement  in  ion  temperature  at  1400  UT — the  peak  of  the 
joule  input.  The  temperature  remains  elevated  even  after  the  joule  heat 
source  terminates  at  1500  UT ,  being  about  150  K  higher  than  the  earlier 
temperatures.  Ttie  temperature  gradually  decreases  and  it  returns  to 
initial  conditions  by  2300  UT . 

The  electron  temperature  shows  the  effect  of  the  changing  Solar- 
Zeni th -Angle  (SZA).  Initially  the  temperature  is  about  2200  K  and  it 
gradually  decreases  with  increasing  SZA  (from  0430  to  0800  UT). 
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FIGURE  5  AVERAGE  ELECTRIC  FIELD  FOR  13  AUGUST  1975.  The  Large  Ex  field 
(~55  mV/m)  at  1400  UT  produced  substantial  joule  heating. 

Tho  ratio  of  the  number  densities  of  0+  ions  to  the  total  number  of 
ions,  [0+]/N  ,  is  used  to  describe  the  ion  composition  vs.  altitude. 
Figure  8  shows  contours  of  L0+|/N^  as  a  function  of  time  for  the  13 
August  event.  The  transition  altitude,  the  altitude  whore  half  the  ions 
are  atomic  (|0*|/N  =  0.5),  is  seen  to  be  near  190  km  at  the  beginning 

of  the  observation.  However,  at  1400  UT ,  the  peak  of  the  heating  event, 
this  altitude  increases  by  50  km — molecular  ions  dominate  below  240  km. 
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FIGURE  6  HEIGHT-INTEGRATED  JOULE  HEATING  FOR  13  AUGUST  1975 


An  important  consequence  of  this  alteration  of  ion  composition  is 
an  immediate  decrease  of  electron  density.  Molecular  ions,  O*  and  N0+ , 
recombine  much  more  rapidly  than  atomic  ions  and  the  result  is  a  decay 
of  the  lower  F  region. 

Electron  density  profiles,  measured  before,  during,  and  after  the 
joule  heating  event  are  shown  in  Figure  9.  The  F-region  depletion  can 
easily  be  seen.  The  critical  frequency  prior  to  the  depletion  was  about 
4.1  MHz.  During  the  heating,  the  critical  frequency  dropped  to  2.7  MHz 
and  30  minutes  later,  it  is  3.5  MHz, 

To  aid  in  the  interpretation  of  these  measurements,  we  performed  a 
numerical  analysis  in  which  we  varied  parameters  such  as  electron  density, 
production  rate,  and  joule  heating  rate.  The  resulting  effects  on  the  ion 
composition  agree  with  the  observed  data,  and  are  described  below. 
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FIGURE  7  ELECTRON  AND  ION  TEMPERATURES  AT  277  km  FOR  13  AUGUST  1975 

B.  Numerical  Modeling  of  Transition  Altitude  Variation 

The  initial  conditions  include  quiet  day  and  night  parameters  in 


order  to  establish  background  composition  profiles.  The  change  in  ion 
composition  can  be  examined  with  the  aid  of  an  expression  derived  using 
the  following  simplified  system  of  reactions  and  continuity  equations: 
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FIGURE  8  ION  COMPOSITION  CONTOURS  FOR  13  AUGUST  1975.  At  a  fixed  altitude 
the  ratio  of  0+  to  total  ionization  decreased  markedly  during  the  joule  heating 
event  and  increased  slowly  thereafter.  The  transition  altitude  computed  during 
the  event  is  a  lower  bound  with  a  5  km  uncertainty. 
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FIGURE  9 


ELECTRON  DENSITY  PROFILES  DURING 
THE  JOULE  HEATING  EVENT  ON  13 
AUGUST  1975.  The  F -region  density 
decreases  during  the  event  and  recovers 
as  the  heating  subsides. 


=  q(0+)  -  Y1[02][0+]  -  V2[N2][0+]  =  0 


^  =  Y2[N2][0+]  +  Y3[N2][0]  -  o2[NO+]Ne  =  0 


a[ot] 


-3—  =  q(02)  +  Y1[02][0  ]  -  c1[02]Ne  =  0 


3[<!  +  +  + 

=  q(N2)  -  Y3[N2][0]  -  a3[N2]Ne  =  0 


where  q(x)  is  the  production  rate  of  species  x.  If  steady  state  is 
assumed,  the  number  densities  of  the  various  ions  are  expressed  as: 


,  Y2[N  ][o+]  y3[n+][o] 

[NcT]  =  2— 2_, -  +  3 
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+.  q(N2) 
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Substitution  of  the  above  into  the  expression  for  charge  neutrality, 


N  =  10+]  +  lNO+ ]  +  10+] 


and  solving  for  [0+]/Nc,  gives 
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and  ?!  =  total  ion  production  rate  [Jones  and  Rees,  1973].  The  effects 
on  the  transition  altitude  of  varying  the  atmospheric  density,  electron 
density,  production  rate,  and  temperature  can  then  be  assessed. 


1  .  Quiet  Conditions — Diurnal  Variation 

The  model  atmosphere  used  in  all  cases  was  taken  from  Banks 
and  Kockarts  for  an  exospheric  temperature  of  1000  K,  and  the  ionospheric 
parameters  were  obtained  from  Cliatanika  radar  data.  The  model  atmospheric 
and  ionospheric  data  are  given  in  Table  3.  The  calculated  composition 
profiles  are  shown  in  Figure  10.  The  difference  in  the  profiles  is  most 
pronounced  above  180  km.  Above  that  altitude,  the  daytime  electron  den¬ 
sity  is  at  least  an  order  of  magnitude  larger  than  the  nighttime  density. 
The  transition  altitude  increases  by  about  20  km  at  night. 
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Tabic  3 

[ODEL  ATMOSPHERE  AND  IONOSPHERIC  DATA  FOR  QUIET  DAY  AND  NIGHT 

Quiet  Day 


Alt 

N 

e 

B 

T 

e 

[°2] 

[0] 

ln2i 

(km) 

(cm  3) 

B 

(K) 

(cm  3) 

(cm  3) 

(cm  3) 

140 

7 .5E4* 

596 

650 

4.40E9 

2.44E10 

5.37E10 

160 

1 .0E5 

765 

1000 

1 . 20E9 

1 .13E10 

1  . 67E10 

180 

1 .2E5 

861 

1600 

4.43E8 

6.45E9 

6.88E9 

200 

1  .  5E5 

916 

1800 

1 .87E8 

4.07E9 

3 . 22E9 

220 

2.0E5 

948 

2000 

8.48E7 

2.69E9 

1 . 60E9 

240 

3.0E5 

967 

2000 

4.00E7 

1 .83E9 

8.26E8 

260 

2.5E5 

979 

2000 

1 .93E7 

1  . 26E9 

4.36E8 

280 

2 . 45E5 

986 

2000 

9.43E6 

8.80E8 

2.33E8 

300 

2.0E5 

991 

2000 

4.67E6 

6.20E8 

1  . 26E8 

Quiet 

Night 

Alt 

N 

c 

B 

LV 

10] 

ln2] 

(km) 

,  -3x 

(cm  ) 

Bj 

,  -3* 

(cm  ) 

(cm  ) 

(cm  3) 

140 

8.6E4 

596 

596 

4 . 40E9 

2.44E10 

5.37E10 

160 

7.0E4 

765 

765 

1 . 20E9 

1  . 13E 10 

1  . 67E10 

180 

5  .6E4 

861 

861 

4.43E8 

6.45E9 

6 . 88E9 

200 

5.0E4 

916 

916 

1 . 87E8 

4.07E9 

3 . 22E9 

220 

3  ,4E4 

948 

943 

8.48E7 

2.69E9 

1  . 60E9 

240 

3.0E4 

965 

965 

4.00E7 

1  . 83E9 

8.26E8 

. 5E4  =  7.5  x  10 
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FIGURE  10  MODEL  COMPOSITION  PROFILES  SHOWING  DIURNAL 
VARIATION 


2  .  Active  Conditions — Particle  Precipitation 

The  effects  of  electron  precipitation  that  are  considered  are 

enhanced  N  ,  electron  temperature,  and  vibrational  temperature  of  . 

The  atmospheric  and  ionospheric  parameters  arc  given  in  Table  4.  The 

first  case  considered  uses  temperature  and  N^  profiles  from  data  taken 

on  18  February  1976.  The  particle  energy  input  rate  was  approximately 
2 

20  ergs/cm  -s .  The  composition  profiles  shown  in  Figure  11  indicate 
that,  relative  to  the  background  profile  (quiet  night),  the  transition 
altitude  is  lower  by  more  than  20  km.  In  fact,  the  transition  altitude 
is  lower  in  this  case  than  it  was  using  daytime  conditions. 

The  effects  of  enhanced  vibrational  temperatures  of  are 
examined  using  the  above  (active  night)  model.  According  to  Sclnink  and 
Banks  [1975],  the  vibrational  temperatures  during  an  auroral  substorm  can 
be  in  the  range  1400  to  2200  K.  These  temperatures  will  increase  the  N 

& 

and  O  reaction-rate  coefficient  (y  )  by  1.5  to  3.6  times  the  ground- 

state  coefficient.  Consequently,  the  composition  profiles  were  recomputed 

applying  those  factors  to  y  .  The  results  arc  shown  in  Figure  12.  There 

is  a  substantial  increase  in  the  transition  altitude  for  the  case  when 

T  =  2200  K. 
v 


Table  4 

MODEL  ATMOSPHERE  AND  IONOSPHERIC  DATA  FOR  ACTIVE  (AURORA)  NIGHT 


Alt 

N 

c 

T 

l 

T 

e 

l°2] 

[0] 

"V 

(km) 

(cm  3 ) 

(K) 

(K) 

,  “3* 

(cm  ) 

(cm  3) 

(cm  3) 

140 

2 . 0E5 

596 

750 

4 . 40E9 

2.44E10 

5 . 3  7E 1  0 

160 

1 .3E5 

765 

1120 

1 . 20E9 

1  . 13E10 

1  . 67E10 

180 

l  .  5E5 

861 

1300 

4 . 43E8 

6.45E9 

6 . 88E9 

200 

1  .  8E5 

916 

i  800 

l  .  87E8 

4.07E9 

3 . 22E9 

220 

1  .  9E5 

948 

2000 

8.48E7 

1  . 83E9 

1 . 60E9 
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FIGURE  11  MODEL  COMPOSITION  PROFILES  SHOWING  THE  EFFECTS 


OF  AURORAL  IONIZATION 


FIGURE  12  COMPOSITION  PROFILES  SHOWING  THE  EFFECT  OF 
ENHANCED  N2  VIBRATIONAL  TEMPERATURE.  The 
model  from  Table  4  was  used  and  the  reaction-rate 
coefficient  for  No  +  0+  was  enhanced. 


Effect  of  Joule  Heat  Input 


The  atmospheric  model  and  ionospheric  data  used  to  evaluate 
the  effect  of  joule  heat  are  given  in  Table  5.  The  ionospheric  data — 
i.e.,  electron  density,  and  electron/ion  temperatures  were  taked  from 
13  August  1975  data.  The  composition  profile  (Figure  13)  corresponding 
to  the  active  data  indicates  that  more  molecular  ions  exist  at  all  alti¬ 
tudes  as  compared  to  the  quiet  period.  The  transition  altitude  increased 
by  20  km.  The  active  period  was  characterized  by  elevated  ion  tempera¬ 
ture  and  decreased  electron  density;  both  factors  contribute  to  the 
increase  of  the  transition  altitude. 

Table  5 

MODEL  ATMOSPHERE  AND  IONOSPHERIC  DATA  USED  TO  EVALUATE 
THE  EFFECT  OF  JOULE  HEAT 
Active  Day  (Joule  Heating) 


Alt 

N 

e 

T. 

l 

T 

c 

l°2J 

10  J 

IV 

(km) 

,  -3. 

(cm  ) 

(K) 

(K) 

(cm  3) 

(cm  3) 

,  -3v 

(cm  ) 

160 

2.0E5 

830 

965 

1 .20E9 

1.13E10 

1  .67E10 

190 

1  .  2E5 

1050 

1400 

2.85E8 

5 . 08E9 

4.66E9 

220 

1  .  5E5 

1200 

2000 

8.48E7 

2  .69E9 

1  . 60E9 

276 

1  .OF,  5 

1800 

2000 

9.60E6 

9.00E8 

2.50E8 

Quiet  Day 
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(km) 

N 

c 

,  "3* 

(cm  ) 

■ 

T 

e 

(K) 

l°2] 

(cm  ) 

[0] 

-3 

(cm  ) 

[n2] 

,  ~3s 

(cm  ) 

160 

2.0E5 

765 

965 

1 . 20E9 

1 . 13E10 

1 . 67EI 0 

190 

2.3E5 

1100 

1784 

2.85E8 

5.08E9 

4.66E9 

220 

3.5E5 
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2000 

8.48E7 

2.69E9 

1  . 60E9 
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FIGURE  13  COMPOSITION  PROFILE  CORRESPONDING  TO  THE  JOULE 
HEAT  INPUT  COMPARED  TO  THE  QUIET  PERIOD 
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Summary  of  Numerical  Modeling 

We  have  just  shown  that  the  ion  composition  at  a  given  altitude 
changes  as  a  function  of  the  production  and  loss  rates.  Factors  that 
contribute  to  the  enhancement  of  the  production  rate  tend  to  increase 
the  relative  abundance  of  0+  ions.  These  factors  are  Solar  EUV  radiation 
during  the  day  and  energetic  auroral  electron  bombardment  at  night.  Both 
the  radar  data  and  the  numerical  modeling  indicate  lower  transition  alti¬ 
tudes  during  these  periods. 

Factors  that  contribute  to  the  enhancement  of  the  loss  rate  tend  to 
increase  the  relative  abundance  of  molecular  ions.  These  factors  include 
elevated  ion  and  neutral  temperatures,  and  enhanced  vibrational  tempera¬ 
tures  of  N2  molecules.  The  ion  temperature  and  the  N2  vibrational  tem¬ 
perature  determine  the  rate  coefficient  and  the  neutral  temperature 
defines  the  Ng  number  density. 


IV  CONCLUSIONS 


Impulsive  ionospheric  heating  events  resulting  from  auroral  energy 
input  have  been  observed  by  the  Chatanika  radar.  In  some  cases,  the 
time — and  height — integrated  auroral  energy  deposition  is  comparable  to 
that  expected  from  a  high  altitude  nuclear  event  outside  the  fireball 
region . 

In  the  auroral  event,  virtually  all  ionospheric  parameters  (electron 
density,  ion,  electron  and  neutral  temperatures,  and  ion  and  neutral 
composition)  are  strongly  affected.  This  would  also  be  the  case  for 
extended  regions  outside  a  nuclear  fireball .  The  increased  temperatures 
and  the  resulting  atmospheric  heave  contribute  to  a  change  in  ion  compo¬ 
sition  such  that  molecular  ions  (N0+,  0^)  dominate  over  atomic  ions  (0+) , 

Ct 

Since  molecular  ions  recombine  much  more  rapidly  than  atomic  ions,  the 
net  electron  density  decreases.  The  region  affected  by  this  depletion 
is  larger  than  the  region  affected  by  the  fireball  and  debris. 

The  nuclear  effects  codes  do  not  attempt  to  accurately  model  this 
region.  Consequently,  they  do  not  accurately  predict  ionospheric  effects 
of  nuclear  explosions  beyond  the  fireball  region.  In  fact,  the  ambient 
atmosphere  and  ionosphere  models  are  not  realistic  compared  to  observed 
parameters . 

The  Chatanika  radar  provides  relevant  information  concerning  both 
the  ambient  and  aurorally  disturbed  ionosphere.  These  results  can  be 
used  to  improve  the  nuclear  effects  codes  so  that  the  much  larger  region 
beyond  the  fireball  is  more  realistically  modeled. 
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DOST 

ATTN 

DELHD-N-P 

3  cy  ATTN 

RAAE 

4  cy  ATTN 

TITL 

U.S.  Army  Comm. -Elec.  Engrg.  Instal .  Agency 

ATTN 

CCC-CED-CCO,  W.  Neuendorf 

Defense  Technical  Information  Center 

ATTN 

CCC-EMEO,  W.  Nair 

12  cy  ATTN 

DD 

ATTN 

CCC-EMEO-PED,  G.  Lane 

Field  Coittnand 

U.S.  Army  Cormiunications  Command 

Defense  Nuclear  Agency 

ATTN 

CC-OPS-WR,  H.  Wilson 

ATTN 

FCPR 

ATTN 

CC-OPS-W 

Field  Command 

U.S.  Army  Communications  R&D  Command 

Defense  Nuclear  Agency 

Livermore  Division 

ATTN 

DRDCO-COM-RY ,  W.  Kesselman 

ATTN 

FCPRL 

U.S.  Army  Foreign  Science  &  Tech.  Ctr. 

ATTN 

DRXST-SD 

Interservice  Nuclear  Weapons  School 

ATTN 

TTV 

U.S.  Army  Materiel  Dev.  &  Readiness  Cmd. 

ATTN 

DRCLDC,  J.  Bender 

Joint  Chiefs  of  Staff 

ATTN 

J-3,  WWMCCS  Evaluation  Office 

U.S.  Army  Nuclear  &  Chemical  Agency 

ATTN 

C3S 

ATTN 

Library 

ATTN 

J- 37 
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DEPARTMENT  OF  THE  ARMY  (Continued) 

U.S.  Army  Satellite  Conn.  Agency 
ATTN:  Document  Control 

U.S.  Army  TRADOC  Systems  Analysis  Activity 
ATTN:  ATAA-TCC,  F.  Payan,  Jr. 

ATTN:  ATAA-PL 
ATTN:  ATAA-TDC 

DEPARTMENT  OF  THE  NAVY 

Joint  Cruise  Missile  Project  Office 

Department  of  the  Navy 
ATTN:  JCM-G-70 

Naval  Air  Development  Center 

ATTN:  Code  6091,  M.  Setz 

Naval  Air  Systems  Command 
ATTN:  PMA  271 

Naval  Electronic  Svstems  Command 
ATTN:  PME  117-20 
ATTN:  PME  117-211,  B.  Kruger 
ATTN:  PME  117-2013,  G.  Burnhart 
ATTN:  PME  106-13,  T.  Griffin 
ATTN:  Code  3101,  T.  Hughes 
ATTN:  PME  106-4,  S.  Kearney 
ATTN:  Code  501 A 

Naval  Intelligence  Support  Ct. . 

ATTN:  NISC-50 

Naval  Ocean  Systems  Center 

ATT N:  Code  532,  J.  8icke) 

ATTN:  Code  5322,  M.  Paulson 
3  cy  ATTN:  Code  5324,  W.  Moler 

Naval  Research  Laboratory 

ATTN:  Code  7550,  J.  Davis 
ATTN:  Code  6700,  T.  Coffey 
ATTN:  Code  6780,  S.  Ossakow 

ATTN:  Code  7500,  B.  Wald 

Naval  Space  Surveillance  System 
ATTN:  J.  Burton 

Naval  Surface  Weapons  Center 
ATTN:  F-31 

Naval  Surface  Weapons  Center 

ATTN:  Code  F-14,  R.  Butler 

Naval  Telecommunications  Command 
ATTN:  Code  341 

Office  of  Naval  Research 
ATTN:  Code  421 

ATTN:  Code  420 

Office  of  the  Chief  of  Naval  Operations 
ATTN:  OP  604C 
ATTN:  OP  941D 
ATTN:  OP  981N 

Strategic  Systems  Project  Office 

Department  of  the  Navy 
ATTN:  NSP-Z141 
ATTN:  NSP-2722,  F.  Wimberly 
ATTN:  NSP-43 


DEPARTMENT  OF  THE  AIR  FORCE 

Aerospace  Defense  Command 
Department  of  the  Air  Force 
ATTN:  DC,  T.  Long 

Air  Force  Avionics  Laboratory 
ATTN:  AAD,  W.  Hunt 
ATTN:  AAD,  A.  Johnson 

Air  Force  Geophysics  Laboratory 
ATTN:  OPR-1 ,  J.  Ulwick 
ATTN:  LKB,  K.  Champion 
ATTN:  OPR,  A.  Stair 
ATTN:  PHI,  J.  Buchau 
ATTN:  PHP,  J.  Aarons 
ATTN:  PHP,  J.  Mullen 

Air  Force  Weapons  Laboratory,  AFSC 
ATTN:  SUL 
ATTN:  DVC 

Air  Logistics  Comnand 
Deoartment  of  the  Air  Force 

ATTN:  00-ALC/MM,  R.  Blackburn 

Assistant  Chief  of  Staff 
Intel  1 igence 

Department  of  the  Air  Force 
ATTN:  INED 

Assistant  Chief  of  Staff 
Studies  &  Analyses 
Department  of  the  Air  Force 
ATTN:  AF/SASC ,  G.  Zank 
ATTN:  AF/SASC,  W.  Adams 

Ballistic  Missile  Office 
Air  Force  Systems  Command 
ATTN:  MNNH,  M.  Baran 
ATTN:  MNNL ,  S.  Kennedy 
ATTN:  MNNH 

Deputy  Chief  of  Staff 
Operations  Plans  and  Readiness 
Department  of  the  Air  Force 
ATTN:  AFXOKCD 
ATTN:  AFXOXFD 
ATTN:  AFXOKS 
ATTN:  AFXOKT 

Deputy  Chief  of  Staff 
Research,  Development,  &  Acq. 
Department  of  the  Air  Force 
ATTN:  AFROS 
ATTN:  AFRDSP 
ATTN:  AFRDQ 
ATTN:  AFRDSS 

Electronic  Systems  Division 
Department  of  the  Air  Force 
ATTN:  DCKC,  J.  Clark 

Electronic  Systems  Division 
Department  of  the  Air  Force 
ATTN:  XRW ,  J.  Deas 

Electronic  Systems  Division 
Department  of  the  Air  Force 
ATTN:  YSEA 

ATTN:  YSM,  J.  kobelski 
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DEPARTMENT 

OF  THE  AIR  FORCE  (Continued) 

OTHER  GOVERNMENT  AGENCIES  (Continued) 

Foreign  Technology  Division 

Department  of  Commerce 

Air  Force  Systems  Command 

National  Bureau  of  Standards 

ATTN 

SDEC,  A.  Oakes 

ATTN 

Security  Officer  for  R.  Moore 

ATTN 

NIIS  Library 

ATTN 

TQTD,  B.  Ballard 

Department  of  Commerce 

National  Oceanic  &  Atmospheric  Admin. 

Headquarters  Space  Division 

Environmental  Research  Laboratories 

Air  Force  Systems  Conmand 

ATTN 

R.  Grubb 

ATTN 

SZJ,  L.  Doan 

ATTN 

SZJ,  W.  Mercer 

Institute  for  Telecommunications  Sciences 

National  Telecommunications  &  Info.  Admin. 

Rome  Air  Development  Center 

ATTN 

D.  Crombie 

Air  Force  Systems  Command 

ATTN 

L.  Berry 

ATTN 

TSLD 

ATTN 

W.  Utlaut 

ATTN 

OCS,  V.  Coyne 

U.S.  Coast 

Guard 

Rome  Air  Development  Center 

Department  of  Transportation 

Air  Force  Systems  Command 

ATTN 

G-D0E-3/TP54 ,  B.  Romine 

ATTN 

EEP 

DEPARTMENT 

OF  DEFENSE  CONTRACTORS 

Strategic  Air  Command 

Department  of  the  Air  Force 

Aerospace  Corp. 

ATTN 

OOKSN 

ATTN 

T.  Salmi 

ATTN 

XPFS 

ATTN 

N.  Stockwell 

ATTN 

DCXF 

ATTN 

R.  Slaughter 

ATTN 

NRT 

ATTN 

V.  Josephson 

ATTN 

DCX 

ATTN 

F.  Morse 

ATTN 

DCXT,  T.  Jorgensen 

ATTN 

S.  Bower 

ATTN 

DCXT 

ATTN 

D.  Olsen 

ATTN 

I.  Garfunkel 

Headquarters  Space  Division 

Air  Force  Systems  Command 

University  of  Alaska 

ATTN 

SKA,  M.  Clavin 

ATTN 

T.  Davis 

ATTN 

SKA,  C.  Rightmyer 

Analytical 

Systems  Engineering  Corp. 

DEPARTMENT 

OF  ENERGY  CONTRACTORS 

ATTN 

Radio  Sciences 

EG&G,  Inc. 

Analytical 

Systems  Engineering  Corp. 

Los  Alamos 

Division 

ATTN 

Securi ty 

ATTN 

Document  Control  for  D.  Wright 

ATTN 

Document  Control  for  J.  Colvin 

Barry  Research  Communications 

ATTN 

J.  McLaughlin 

Lawrence  Livermore  Laboratory 

ATTN 

Document  Control  for  Technical 

BDM  Corp. 

Information  Dept.  Library 

ATTN 

T.  Neighbors 

ATTN 

L.  Jacobs 

Los  Alamos 

Scientific  Laboratory 

ATTN 

Document  Control  for  P.  Keaton 

Berkeley  Research  Associates,  Inc. 

ATTN 

Document  Control  for  D.  Westervelt 

ATTN 

J.  Workman 

ATTN 

Document  Control  for  R.  Taschek 

Boeing  Co. 

Sandia  Laboratories 

ATTN 

S.  Tashird 

ATTN 

Document  Control  for  Org.  1Z50 

ATTN 

G.  Hall 

W.  Brown 

ATTN 

D.  Murray 

ATTN 

Document  Control  for  3141 

ATTN 

J.  Kenney 

ATTN 

Document  Control  for  Space  Project  Div. 

ATTN 

Document  Control  for  D.  Thornbrough 

University  of  California  at  San  Diego 

ATTN 

Document  Control  for  D.  Dahlgren 

ATTN 

H.  Booker 

Sandia  Laboratories 
Livermore  Laboratory 

ATTN:  Document  Control 
ATTN:  Document  Control 

OTHER  r„)VERNME_NJ  AGENCIES 

Central  Intel  1 igence  Agency 
ATTN:  OSI/PSTD 


for  T.  Cook 
for  E.  Murphey 


Charles  Stark  Draper  Lab.,  Inc. 
ATTN:  D.  Cox 
ATTN:  J.  Gilmore 

Computer  Sciences  Corp. 

ATTN:  H.  Blank 

Comsat  labs 

ATTN:  R.  Taur 
ATTN:  G.  Hyde 
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DEPARTMENT 

OF  DEFENSE  CONTRACTORS  (Continued) 

DEPARTMENT 

OF  DEFENSE  CONTRACTORS 

Cornell  University 

Johns  Hopkins  University 

ATTN 

D.  Farley,  Jr. 

ATTN 

J.  Newland 

ATTN 

T.  Evans 

Electrospace  Systems,  Inc. 

ATTN 

B.  Wise 

ATTN 

H.  Logston 

ATTN 

P.  Komiske 

ATTN 

Document  Librarian 

ESL,  Inc. 

ATTN 

T .  Potemra 

ATTN 

J.  Marshall 

ATTN 

J.  Roberts 

Linkabit  Corp. 

ATTN 

C.  Prettie 

ATTN 

I.  Jacobs 

Ford  Aerospace  &  Communications 

Corp. 

Litton  Systems,  Inc. 

ATTN 

J.  Mattingley 

ATTN 

R.  Grasty 

General  Electric  Co. 

Lockheed  Missiles  &  Space  Co.,  Inc 

ATTN 

M.  Bortner 

ATTN 

W.  Imhof 

ATTN 

R.  Johnson 

General  Electric  Co. 

ATTN 

M.  Walt 

ATTN 

A.  Steinmayer 

ATTN 

C.  Zierdt 

Lockheed  Missiles  &  Space  Co.,  Inc 

ATTN 

S.  Lipson 

ATTN 

Dept  60-12 

ATTN 

D.  Churchill 

General  Electric  Co. 

ATTN 

F.  Reibert 

M. I ,T .  Li ncol n  Lab 

ATTN 

L.  Loughlin 

General  Electric  Company-TEMPO 

ATTN 

D.  Towle 

ATTN 

M.  Stanton 

ATTN 

T.  Stevens 

McDonnell  Douglas  Corp. 

ATTN 

DAS  I  AC 

ATTN 

J.  Moule 

ATTN 

W.  Knapp 

ATTN 

W.  Olson 

ATTN 

D.  Chandler 

ATTN 

N.  Harris 

ATTN 

G.  Mroz 

General  Electric  Tech.  Services 

Co., 

Inc. 

ATTN 

G.  Millman 

Mission  Research  Corp. 

ATTN 

R.  Hendrick 

General  Research  Corp. 

ATTN 

S.  Gutsche 

ATTN 

J.  Ise,  Jr. 

ATTN 

R.  Bogusch 

ATTN 

J.  Garbarino 

ATTN 

F.  Fajen 

ATTN 

D.  Sowle 

GTE  Sylvania,  Inc. 

ATTN 

M.  Cross 

Mitre  Corp 

ATTN 

B.  Adams 

HSS,  Inc. 

ATTN 

A.  Kynriel 

ATTN 

D.  Hansen 

ATTN 

C.  Callahan 

ATTN 

G.  Harding 

IBM  Corp. 

ATIN 

F.  Ricci 

Mitre  Corp 

ATTN 

W.  Hall 

Uni versi ty  of  Ill inoi s 

ATTN 

W.  Foster 

ATTN 

Security  Supervisor  for 

K.  Yeh 

ATTN 

M.  Horrocks 

Institute  for  Defense  Analyses 

Pacific-Sierra  Research  Corp. 

ATTN 

E.  Bauer 

ATTN 

E.  Field,  Jr. 

ATTN 

J.  Aein 

ATTN 

H.  Wolfhard 

Photometries,  Inc. 

ATTN 

J.  Bengston 

ATTN 

I.  Kofsky 

International  Tel.  8  Telegraph  Corp. 

Physical  Dynamics,  Inc. 

ATTN 

G.  Wetmore 

ATTN 

E.  Fremouw 

ATTN 

Technical  Library 

R  &  D  Associates 

JAYCOR 

ATTN 

B.  Yoon 

ATTN 

S.  Goldman 

ATTN 

L.  Delaney 

JAYCOR 

Rand  Corp. 

ATTN 

D.  Carlos 

ATTN 

C.  Crain 

ATTN:  E.  Bedrozian 


Kaman  Sciences  Corp. 

ATTN:  T,  Meagher 
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DEPARTMENT 

OF  DEFENSE  CONTRACTORS  (Continued! 

DEPARTMENT 

OF  DEFENSE  CONTRACTORS 

R  &  D  Associates 

SRI  International 

ATTN 

W.  Karzas 

ATTN 

W.  Jaye 

ATTN 

M.  Gantsweg 

ATTN 

A.  Burns 

ATTN 

R.  Turco 

ATTN 

G.  Smith 

ATTN 

C.  Greifinger 

ATTN 

G.  Price 

ATTN 

H.  Ory 

ATTN 

M.  Baron 

ATTN 

W.  Wright,  Jr. 

ATTN 

D.  Neilson 

ATTN 

C.  MacDonald 

ATTN 

C.  Rino 

ATTN 

B.  Gabbard 

ATTN 

R.  Leadabrand 

ATTN 

R.  Lelevier 

ATTN 

R.  Livingston 

ATTN 

F.  Gilmore 

ATTN 

W.  Chesnut 

3  cy  ATTN 

J.  Kelly 

Riverside  Research  Institute 

ATTN 

V.  Trapani 

Teledyne  Brown  Engineering 

ATTN 

R.  Deliberis 

Rockwell  International  Corp. 

ATTN 

J.  Kristof 

Tri-Com,  Inc. 

ATTN 

D.  Murray 

Santa  Fe  Corp. 

ATTN 

E.  Ortlieb 

TRW  Defense  X  Space  Sys.  Group 

ATTN 

S.  Altschuler 

Science  Appl ications ,  Inc. 

ATTN 

R.  Plebuch 

ATTN 

L.  Linson 

ATTN 

D.  Dee 

ATTN 

C.  Smith 

ATTN 

D.  Hamlin 

Utah  State 

University 

ATTN 

J.  McDougall 

ATTN 

K.  Baker 

ATTN 

D.  Sachs 

ATTN 

L.  Jensen 

ATTN 

E.  Straker 

Visidyne,  Inc. 

Science  Appl ications,  Inc. 

ATTN 

J.  Carpenter 

ATTN 

D.  Divis 

Science  Applications,  Inc. 

ATTN:  SZ 
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